Pseudoscalar perturbations and polarization of the cosmic microwave background 
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We show that models of new particle physics containing massless pseudoscalar fields super-weakly 
coupled to photons can be very efficiently probed with CMB polarization anisotropies. The stochas- 
tic pseudoscalar fluctuations generated during inflation provide a mechanism for converting i5-mode 
polarization to _B-mode during photon propagation from the surface of last scattering. The efficiency 
of this conversion process is controlled by the dimensionless ratio H/{2Tifa), where H is the Hubble 
scale during inflation, and fa^ is the strength of the pseudoscalar coupling to photons. The current 
observational limits on the -B-mode constrain this ratio to be less than 0.07, which in many models 
of inflation translates to a sensitivity to values of fa in excess of 10^* GeV, surpassing the sensitivity 
of other tests. 



Introduction: — Within the last decade, precision ob- 
servations of fluctuations in the cosmic microwave back- 
ground (CMB) have provided us with a powerful probe 
of cosmology, and been the driving force in pinning down 
the parameters in the concordance ACDM model. In re- 
cent times, gains in sensitivity have allowed us to comple- 
ment observations of the temperature fluctuations with 
those in various components of polarization. Much of the 
focus has been aimed at the detection of the so-called 
-B-mode of polarization, as it is an important probe of 
primordial gravitational waves, and thus of the possi- 
ble mechanism for inflation. However, precision probes 
of polarization in various modes also present us with a 
powerful set of tools to test for new degrees of freedom 
with which the CMB photons may interact. In this Let- 
ter, we show that existing upper bounds on the B-mode 
of polarization already place stringent constraints on any 
new light pseudoscalar degrees of freedom through their 
stochastic fluctuations generated during inflation. 

The interactions of a massless (or nearly massless) 
pseudoscalar field a with photons may be parametrized 
in terms of a dimensionful coupling constant /^f^. 



where is the electromagnetic field strength tensor, 
and Ffj_i, its dual. We note that the absence of a potential, 
V{a) = 0, is protected by the shift symmetry of this 
Lagrangian, a a + const, and while it is also natural 
to expect that the pseudoscalar a would have derivative 
couplings to the spins of other particles, this does not 
affect the physics of the CMB. The presence of the photon 
coupling implies that a space-time variation of the field 
a will induce a rotation Atp of the polarization of an 
electromagnetic wave. Provided the photon wavelength 
is much smaller than the characteristic scale on which 
the field a varies , the rotation angle is independent of 



frequency: 



Alp = —Aa. 

J a 



(2) 



In the past, CMB polarization has been used to probe an 
induced rotation angle due to a spatially homogeneous 
pseudoscalar profile a — a{t) 0, with the latest data 
Ul, 01 giving comparable sensitivity to the strongest con- 
straints from polarized extra-galactic radio sources 0]. 
However, the time evolution of the background effectively 
amounts to a bulk violation of parity and/or Lorentz in- 
variance. Here, we consider a scenario perfectly in keep- 
ing with standard cosmology, and focus on a stochastic 
background for the pseudoscalar as would arise from a 
period of inflation. We exploit the fact that all mass- 
less scalars have a universal amplitude for fluctuations 
generated during inflation, 5a — H/{2n), where H is the 
inflationary Hubble parameter at horizon crossing. With- 
out loss of generality, we are allowed to choose Onow = 0, 
and then the additional rotation ip ([2]) for a photon arriv- 
ing from the direction n is simply given by the randomly 
fluctuating value of A{n,T) = a{n,T)/fa at the surface 
of last scattering (LSS), t ~ tlsSi where t is the con- 
formal time. To quantify the effect, we introduce the 
dimensionless parameter Ca'. 
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It is then apparent that any 0(1) probe of these fluctua- 
tions will provide impressive sensitivity to fa on the or- 
der of the inflationary Hubble scale. For a scale-invariant 
source of inflationary perturbations in the field a, the evo- 
lution of the rotation angle ip with the scale factor R{t) 
can be pictured as an effective random walk where each 
efold of expansion corresponds to a separate step, so that 
the cumulative effect can be further enhanced relative to 

©: (A^)2 ~Caln(i?now/i?). 

In what follows, we analyze the generation of the B- 
mode of CMB polarization from the inflationary fluctua- 
tions of a massless pseudoscalar. We use the B-mode 
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spectrum to deduce a precise limit on Ca, and finally 
describe its implications for models of particle physics 
incorporating new pseudoscalar fields. We should note 
that the impact of inflationary axion perturbations has 
previously been discussed in the context of isocurvature 
perturbations on the CMB temperature spectrum Q, 
but this is not directly relevant for non-axionic mass- 
less pseudoscalars. The technical aspects of our analysis 
do nonetheless have parallels with studies of the impact 
on the CMB spectra of stochastic magnetic fields 7], of 
scalar moduli [8[, and of weak lensing, in that the con- 
volution of an additional fluctuating source with the un- 
derlying density perturbations tends to mix the various 
power spectra for temperature and E and i?-mode po- 
larization. 

Generating B-mode polarization: — At the LSS the un- 
derlying spectrum of scalar (and tensor) perturbations 
leaves its imprint on the CMB through anisotropics in 
temperature and linear polarization, the latter generated 
by Thomson scattering. Linear polarization is conve- 
niently described in terms of the Q and U Stokes param- 
eters, and while a single scalar perturbation mode with 
k\\z does not generate f7, it generates Q. We general- 
ize the approach of Zaldarriaga and Seljak [9] where, for 
a mode of momentum k directed along z, the observed 
value of (5 in a given direction fi can be written as an 
integral over conformal time r along the line of sight, 



Q{k, h) = ^(1 - ^2) p e"^5(T)n(fc, r). 



(4) 



Here x — k{To — t) and = k ■ n. The source and vis- 
ibility functions, n(fc,T) and .g(r) reflect the details 
of the generation of polarization via Thomson scatter- 
ing at the LSS. In the presence of a small stochastic 
rotation ([2]), a ?7 component would also be generated, 
U ~ 2Ai/jQ = 2AQ. This rotation also depletes Q, but 
this is a less significant effect. It then follows that for the 
combination of a single Fourier mode of the underlying 
scalar perturbation, with k\\z, and one Fourier mode of 
the pseudoscalar perturbation with arbitrary direction q, 
the result for U can be exressed as 

U{k,q,n) = ^(1 -m') P e^'^+'^y''g{T)nik,T)AA{r,q). 

. (5) 

In this expression, y ~ q{To — t) and v = q ■ n. 
Aa(9,t) is the transfer function, that has a generic 
form for any massless (pseudo)scalar field: AA{q,T) is 
normalised to one on large scales, and oscillates with 
a decaying envelope determined by the matter content 
on sub-horizon scales. At the next step, we calculate 
basis-independent expansion coefficients by performing 
the projection onto spin 2-weighted spherical harmonics, 
asim = - J d^{Y2im + ^-2 ;m)^(")/2, as described in 
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dn^d^kd^qY*^l„^{h) 



X / dT{m^ - (1 + dlfx^)e'''>'+'y''F{T, fc, q). (6) 

"'0 

Since the aBim^s are basis-independent, we are now al- 
lowed to generalize the scalar-pseudoscalar source func- 
tion to modes of arbitrary direction k and q, which are in 
turn determined by the primordial scalar ^(fc) and pseu- 
doscalar S,a{q) perturbations: 

F(t, k, q) = 5(T)n(fc, r)A^(g, T)ak)U{q)- (7) 

The standard assumption that these perturbations are 
gaussian random variables implies 

(r(fci)^(fc2)) = P0(fci)5(=^)(fci-fc2), 

{CAmUm) ^ PA{qi)S^'Hqi-q2), (8) 

where is the (now rather well-measured) primordial 
power spectrum for scalar perturbations. The scalar pri- 
mordial power spectrum sources the CMB anisotropy 
and indeed structure formation, while the pseudoscalar 
power spectrum Pa is, in our approach, directly linked 
to Ca in ^ by the standard inflationary prediction, 



(9) 



Here is the spectral index of the pseudoscalar primor- 
dial power spectrum, with tIq = 1 being scale invariant. 
Note that the factor of Air reflects the normalization con- 
vention for Fourier transforms in Q that we follow in this 
paper. 

To proceed, we expand the exponential e"'^+*^'^ in ([6]) 
into a basis of products of spherical harmonics, and using 
their orthogonality and completeness, perform all the an- 
gular integrals. Deferring the remaining details J^, we 
obtain 
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X {2h + l){2h + l){ I 0^ 



m , / 1 , / 2 



with the generalized transfer function, 



z,„(fc,g)|2, (10) 



3 f^" 

Ai^i^rn(k,q) = drg{T)ji^{x)ji^{y) 

|l±|j|ij - m^^ AA{r, g)n(r, k). (11) 

The sum over m can be straightforwardly computed, 
while the summation range for li and I2 is restricted 
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FIG. 1: We display the Ci's for the E-mode (green - dot-dash) 
and induced B-mode polarization (red - solid). The fiducial 
value Ca = 4.2 x 10~^ was chosen to match the current upper 
bound on the B-mode using data from QUaD. 95% upper 
limits from QUaD (stars) and 2a upper limits from WMAP5 
(polygons) are shown for comparison. We also exhibit the 
Cbi's induced by lensing of the i5-mode (cyan - long-dash), 
and the contribution from primordial gravitational waves with 
r = 0.14 (blue - dash), corresponding to the choice H14 = 1. 



to 1^1 — ^2 1 < I < ^1 + 1-2, which is enforced by the 
Wigner 3j-symbol. It is important to note that the cross- 
correlations TB and EB vanish on account of the overall 
conservation of parity. 

At this point we assume a spectrum of scalar fluctua- 
tions with scalar spectral index n = 0.963 chosen accord- 
ing to the best fit model for the WMAP5 CMB data, and 
take the pseudoscalar spectral index Ua to have the same 
value. We used our own code based on CMBf ast [ij] to 
numerically calculate the functions n(fc,T) and A^(t, g) 
for the WMAP5 best fit model, and then to compute Cbi 
using Eqs. (P)- pT|) . Fig. [1] displays the results for Cbi 
for a fiducial choice of Ca = 4.2 x 10~^. For comparison, 
the plot also shows Cei as well as the B-modes gener- 
ated by primordial tensor perturbations with r = 0.14 
and by lensing of the i?-mode. The qualitative form of 
Cbi induced by pseudoscalar perturbations can be un- 
derstood by looking at the dominant regions in the li 
and I2 sums. For large I, we find that li ^ I, while the 
sum over I2 is effectively truncated at a lower value of 
order l!^'^^ ~ {tq — tlss)/tlss 50 (higher values of I2 
contribute no more than about 2% to the Cbi's on small 
scales). Thus for large / ~ 1000 the induced B-mode 
closely tracks the underlying E-mode. For lower values 



of I 0(1), both li and I2 saturate at higher scales so 
the B-mode is somewhat larger. Finally, the overall scale 
of the oscillations in Cbi is slightly suppressed as is to be 
expected from convoluting the underlying i?-mode source 
with a gaussian random variable. 

It is important to keep in mind that our approximation, 
which assumes a small rotation angle, may break down 
once Cbi becomes comparable to Cei- Consequently, for 
setting constraints on Cbi we choose I in the interval 
100-1000, where the most recent experimental results of 
QUaD [U] probe the B-modes well below the detected 
E-mode level, Cbi O(O.l) x Cei- The QUaD limits, 
shown in Fig. [U impose a stringent constraint on c^: 

Ca < 4.2 X 10"^ =^ fa> 2.4 X lO^'' GeV x Hu, 

(12) 

where we have introduced H14 = H/IQ^^ GeV, a normal- 
ization inspired by the fact that H14 ^ 0(1) is close to 
the maximal inflationary value of H allowed by observa- 
tions. The inflationary Hubble scale can be traded for the 
tensor-to-scalar ratio r, commonly used to parametrize 
the strength of all massless perturbations including grav- 
itational waves, r — 0.14 x Hf^. Note that from Eq. ([3]), 
the conclusion that Ca ^ 1 justifies a posteriori our per- 
turbative treatment of the rotation of polarization. The 
constraint (|12p is the main result of this Letter, and in 
what follows we discuss its implications. 

Implications for particle physics: Given the conven- 
tional picture of infiationary cosmology, the constraint 
obtained above applies to massless pseudoscalars (or al- 
most massless with a mass below the Hubble scale at 
decoupling) and it is important to consider how such 
new low-energy degrees of freedom may naturally arise. 
Recall that massless pseudoscalar fields coupled to the 
operator G^vG^^" in QCD, namely axions, resolve the 
strong CP problem in a natural way [l3| while gaining 
an anomaly-induced mass ma ''^ m-nfTr/ fa- While this 
mass is large compared to the scales relevant here, this 
mechanism for inducing a mass is unique, and thus should 
two or more pseudoscalars couple to the QCD anornaly, 
only one linear combination would become massive [14| . 
Schematically, below the QCD scale, such models lead to 
an effective Lagrangian of the form. 
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Qi a2_ 
2/1 2f2, 

CQCBa + T^F^^F'^-', (13) 



where besides the Lagrangian /Cqcdo for the QCD-axion, 
one has a massless pseudoscalar a as part of C-ya in ([1]) 

with fa^ = (.92 7/2 - gi/fi)/ ^/gl+gl- in the absence of 
any special reasons for cancellation, the generic case is 
f~^ ^ 0. Although we refrain from assessing the likeli- 
hood of new high-energy physics leading to two or more 
light pseudoscalar fields, we can refer to various scenar- 
ios in string theory where multiple pseudoscalar moduli 
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are a generic prediction (see e.g. [l5|)- Therefore, the 
Lagrangian ([1]) can easily originate from a more funda- 
mental theory with multiple degrees of freedom. Note 
that it would take a significant increase in sensitivity be- 
fore the limit in Eq. could start probing fa at the 
scales most relevant for string theory, fa ^ 10^^ GeVi?i4. 

Existing constraints on the coupling fa arise via several 
known mechanisms. Most notably, the CAST experiment 
directly limits the emission of pseudoscalars from the so- 
lar interior and has obtained the limit fa > 2x 10^° GeV 
[l6| . which is quite competitive with stellar constraints 
that typically require fa 10^^ GeV [17|. A massless 
pseudoscalar-photon coupling in the intergalactic mag- 
netic field provides an additional means of probing fa 
[lit which is, however, strongly dependent on the as- 
sumptions concerning the strength of magnetic field, its 
redshift behavior and the number density of free elec- 
trons. Thus we observe that for H14 ^ 0(1), the con- 
straint (jl2|) derived in this paper is more stringent than 
any other constraint on massless pseudoscalars by at least 
two orders of magnitude. 

Concluding Remarks: — Despite all the evidence for 
physics beyond the Standard Model in the neutrino sec- 
tor and in cosmology, there is no compelling need for new 
degrees of freedom at or below the electroweak scale. In- 
deed, neutrino masses and dark matter can arise from 
ultra-heavy degrees of freedom, while dark energy is cur- 
rently consistent with being just a cosmological constant. 
Nonetheless, many models of new high-energy physics 
do lead to new light degrees of freedom whose mass 
may be protected by symmetry, as is the case for the 
pseudoscalars considered here, and we have shown that 
CMB physics can potentially be a very important probe 
of such models. The stochastic fluctuations of a light 
pseudoscalar field generated during inflation can induce 
the conversion of i?-mode to S-mode polarization in the 
CMB, and with the current upper bounds this already 
leads to stringent constraints ([H]). It is anticipated that 
forthcoming experiments such as BICEP will significantly 
improve the sensitivity to H/{2Tifa) provided the back- 
ground due to lensing can be handled appropriately, and 
thus these experiments may provide the primary sensi- 
tivity to new physics of this type. 

The search for primordial gravitational waves using B- 
modes has intrinsic limitations in that the effect scales 
as (H/Mpi)^, and so for H below 10" GeV the de- 
tection of tensor modes in the CMB becomes problem- 
atic. In this respect, the search for pseudoscalar fluctu- 
ations provides an independent motivation for studying 
B-modes, which is less predicated on the scale of infla- 
tion being large. Indeed, with the existing constraints 
on fa [16], _B-modes may be generated at an observ- 
able level within a wide class of inflationary models with 
10" GeV < H < IQi'* GeV. However, since both pseu- 
doscalars and gravitational waves (when combined with 
the background contribution from lensing) can have sim- 



ilar power spectra up to normalization, it will be impor- 
tant to investigate whether the two can be distinguished 
observationally [l3|- 
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